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Abstract 

The  flow  distribution  in  fuel  cells  has  an  important  influence  on  both  the  power  density  and  efficiency  of  fuel  cell  systems.  In  order 
to  effectively  utilize  the  area,  flow  distribution  should  be  as  homogeneous  as  possible.  In  addition,  pressure  losses  should  be  minimized 
with  regard  to  the  power  demand  of  auxiliary  components  as  pumps  and  compressors.  In  polymer  electrolyte  fuel  cells  (PEFCs)  and  direct 
methanol  fuel  cells  (DMFCs)  the  flow  field  is  in  direct  contact  with  the  diffusion  layer.  The  main  task  of  the  diffusion  layer  is  to  distribute 
the  reactants  from  the  flow  field  towards  the  catalyst  layer.  To  prevent  diffusion  overvoltages,  the  diffusion  layer  is  in  general  highly  porous 
and  provides  high  fluxes  of  the  reactants.  Consequently,  the  flow  distribution  in  the  flow  field  can  be  superpositioned  by  a  flow  in  the 
diffusion  layer.  In  this  paper,  we  discuss  the  interaction  between  the  diffusion  layer  and  the  flow  field.  Experimentally,  we  characterized 
different  diffusion  layers  with  regard  to  their  diffusion  properties  as  well  as  different  flow  fields.  Additional  simulation  studies  help  to 
understand  the  processes  and  to  determine  suitable  combinations  of  flow  fields  and  diffusion  layers. 
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Keywords:  Diffusion  layer;  Polymer  electrolyte;  Fuel  cells 


1.  Introduction 

Fuel  cells  are  promising  energy  converters  since  the  effi¬ 
ciency  achievable  is  higher  than  in  power  plants  or  internal 
combustion  engines. 

The  bipolar  plates  in  fuel  cell  stacks  have  to  fulfill  var¬ 
ious  requirements  such  as  homogeneous  flow  distribution, 
current  transport  and  product  water  removal.  Furthermore, 
the  bipolar  plates  are  in  direct  contact  with  the  membrane 
electrode  assemblies  (MEAs)  so  that  they  have  to  be  highly 
conductive  and  chemically  inert.  For  economic  reasons  and 
due  to  market  demands,  they  are  required  to  be  light-weight 
and  to  be  manufactured  at  low  costs. 

In  particular  PEFCs  are  interesting  for  automotive  and 
portable  applications  because  of  their  low  operating  temper¬ 
ature.  In  this  cell  types,  the  cathodes  have  a  great  impact  on 
the  performance  of  the  cell.  To  reduce  the  energy  demand 
of  the  auxiliary  components  such  as  compressors  or  pumps, 
the  flow  rate  and  the  pressure  drop  should  be  as  low  as 
possible.  The  basic  processes  of  both  PEM  and  DMFC  fuel 
cells  are  shown  in  Fig.  1 .  In  addition  to  the  electrochemical 
oxidation  and  reduction  processes  the  performance  of  the 
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DMFC  is  affected  by  the  methanol  permeation  from  the 
anode  to  the  cathode  where  the  methanol  is  chemically  ox¬ 
idized.  Although  PEM  and  DMFC  have  similar  theoretical 
open  circuit  voltages  Eq  (Table  1)  the  methanol  perme¬ 
ation  lowers  the  DMFC  cell  voltage  by  several  hundreds 
of  mV. 

Generally,  a  fuel  cell  stack  consists  of  several  single  cells 
connected  in  series  by  bipolar  plates.  In  the  same  time,  the 
flow  fields  of  the  bipolar  plates  provide  the  MEAs  with 
reactants.  Fig.  2  gives  a  more  detailed  view  inside  a  single 
cell.  Each  flow  field  is  in  direct  contact  with  the  membrane 
electrode  assembly  which  consists  of  the  catalyst  coated 
proton  conducting  membrane  and  the  diffusion  layers.  To 
provide  high  fluxes  of  reactants  towards  the  catalyst  layer 
the  diffusion  layers  are  highly  porous.  In  addition,  product 
water  removal  is  made  by  a  fraction  of  hydrophobic  pores 
which  is  dependent  on  the  PTFE  content  inter  alia.  The 
task  of  the  flow  field  is  to  provide  a  homogeneous  distri¬ 
bution  of  the  reactants  on  the  whole  cell  area.  An  effective 
use  of  the  whole  flow  field  is  important  with  regard  to  the 
power  output,  the  specific  volumetric  power  density  and  at 
least  the  utilization  of  the  catalyst  layer.  DMFC  stacks  in 
the  power  range  from  several  watts  to  the  kW  class  have 
recently  been  reported  [1-3]. 

In  addition  to  the  flow  homogeneity,  the  pressure  drop 
from  the  inlet  manifold  to  the  outlet  should  be  as  low  as 
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Nomenclature 

a 

anode 

^cell 

cell  area  (m2) 

^mem 

active  area  of  the  single  cell  (m2) 

Am 

cross-sectional  area  of  the 

meander  channel  (m2) 

t>M 

width  of  the  meander  channel  (m) 

c 

cathode 

CMeOH 

methanol  concentration  (mol/1) 

t/hyd 

hydraulic  diameter  (m) 

£o 

open  circuit  voltage  (V) 

F 

Faraday  constant  96,485  As/mol 

hu 

height  of  the  meander  channel  (m) 

in 

at  the  inlet 

i 

cell  current  (A) 

k 

permeability  (m2) 

L 

coupling  matrix 

Cm 

length  of  the  meander  (m) 

4Ca;r 

air  mass  28.96  kg/mol 

N 

number  of  parallel  channels 

out 

at  the  outlet 

p 

pressure  (bar) 

p 

vector  containing  the  local  pressures 

r 

radius  (m) 

Rin 

inner  diameter  (m) 

^out 

outer  diameter  (m) 

Re 

Reynold’s  number 

s 

gap  width  (m) 

Um 

circumfence  of  the  meander  channel  (m) 

T’difT 

flow  velocity  in  the  diffusion  layer  (m2) 

Pspec 

specific  volume  flow  (m3/m2s) 

V 

volume  flow  (m3/s) 

X 

molar  fraction  (mol/mol) 

Greek  letters 

X 

fraction 

h 

dynamic  viscosity  (kg/ms) 

X 

stoichiometric  flow  factor 

V 

kinematic  viscosity  =  r)l p  (Pa  s) 

p 

density  (kg/m3) 

i jr 

flow  homogeneity 

possible  because  of  the  power  consumption  of  the  auxiliary 
components  (e.g.  of  the  air  compressor). 

Different  flow  field  structures  are  well  known  for  PEM 
and  DMFC  fuel  cells  (Fig.  3)  [4—8].  In  principle,  the  channel 


structure  offers  a  homogeneous  flow  distribution  with  low 
pressure  drops.  But  on  the  other  hand  the  formation  of  liquid 
water  droplets  at  the  cathode  can  flood  one  or  more  channel 
with  the  consequence  of  stopping  the  gas  flow  there.  In  pre¬ 
vious  measurements  with  100  cm2  laboratory  test  cells  we 
observed  channel  flooding  after  a  certain  time  of  constant 
operation  (approximately  10  min)  which  leads  directly  to  a 
voltage  decrease.  The  decrease  can  be  stopped  only  by  re¬ 
moving  the  water  frequently  from  the  channels  by  purging 
the  cathode  [6].  In  single  cells  this  may  not  be  a  problem, 
but  in  stacks  the  probability  of  channels  being  blocked 
by  water  droplets  is  elevated  simply  due  to  the  increased 
number  of  cells.  This  requires  more  often  compressor 
power  consuming  purging  procedures  affecting  the  overall 
efficiency. 

Therefore,  special  attention  must  be  drawn  to  a  flow  field 
with  no  flooding  effects  and  low  pressure  drops. 

With  regard  to  the  flooding  the  interdigitated  structure  is 
beneficial.  The  flow  is  pressed  from  the  inlet  to  the  outlet 
directly  through  the  diffusion  layer.  This  structure  provides 
also  the  regions  of  the  catalyst  layer  which  are  situated  un¬ 
derneath  the  rib  with  sufficient  amounts  of  reactants.  As  a 
drawback,  the  pressure  drop  from  the  inlet  to  the  outlet  is 
relatively  high  [4,6]. 

The  third  flow  field  (c)  shown  in  Fig.  3  is  the  meander 
structure  where  one  or  more  parallel  meanders  guide  the  flow 
across  the  active  area  of  the  membrane-electrode-assembly. 
Increasing  the  number  of  meander  channels  decreases  the 
pressure  drop  because  the  total  volume  flow  is  distributed 
and  the  length  of  each  meander  decreases  with  their  num¬ 
ber,  but  on  the  other  hand  too  many  meanders  lead  to  the 
same  situation  described  for  the  parallel  channel  flow  field, 
i.e.  partially  meander  channels  can  be  blocked  by  water 
droplets. 

1.1.  Flow  interaction  between  the  diffusion  layer 
and  the  flow  field 

In  both  PEM  and  DMFC  fuel  cells  the  flow  field  is  in 
direct  contact  with  the  diffusion  layer  which  distributes  the 
reactants  from  the  flow  field  towards  the  catalyst  layer.  To 
avoid  diffusion  overvoltages  the  diffusion  layer  is  in  gen¬ 
eral  highly  porous  and  provides  high  fluxes  of  the  reactants. 
Consequently,  the  flow  distribution  in  the  flow  field  can  be 
superpositioned  by  a  flow  in  the  diffusion  layer  which  is 
schematically  indicated  in  Fig.  4.  In  this  paper,  we  focus  on 
the  flow  interaction  mainly  in  meander  structures. 

In  the  following  the  different  approaches  for  our  study 
are  described.  In  experiments  simulating  the  processes  in 


Table  1 

Reactions  and  voltages  of  PEM  and  DMFC 


Type 

Overall  reaction 

—  AH0  (kJ/mol) 

-AG0  (kJ/mol) 

E0  (V) 

rh h  (AG/AH) 

PEM 

Ht  +  0.502  -*■  PEOliquid 

285.8 

237.4 

1.229 

0.829 

DMFC 

1.502  +  CH3OH  -»  C02  +  2H20iiquid 

726.3 

702.4 

1.214 

0.967 

1  2  3  4  5 


Fig.  2.  Membrane  fuel  cell.  1:  anode  bipolar  plate  with  flow  field  for 
electrical  contact  and  distribution  of  reactant;  2:  anode  diffusion  layer;  3: 
catalyst  coated  membrane;  4:  cathode  diffusion  layer;  5:  cathode  bipolar 
plate. 

fuel  cells  we  visualized  the  flow  distribution  in  differ¬ 
ent  combinations  of  meanders  and  diffusion  layers.  In  a 
second  step,  a  model  was  developed  describing  the  flow 
processes.  The  results  of  parameter  studies  serve  to  de¬ 
termine  suitable  combinations  of  flow  fields  and  diffusion 
layers. 


Fig.  4.  Scheme  of  the  interaction  between  the  flow  field  and  the  diffusion 
layer.  The  diffusion  layer  is  in  direct  contact  with  the  flow  field.  Due  to 
the  pressure  gradient  in  the  flow  field,  there  is  an  additional  flow  in  the 
diffusion  layer  which  is  superpositioned  by  a  flow  in  the  diffusion  layer. 
The  resulting  overall  flow  distribution  is  the  sum  of  both  mass  flows. 

2.  Characteristic  parameters  of  flow  distribution 

This  section  briefly  describes  definitions  and  characteris¬ 
tic  parameters  of  the  flow  in  fuel  cells. 

2.1.  Flow  homogeneity 

Cell  areas  supplied  with  low  amounts  of  reactants  are  af¬ 
fected  by  diffusion  overvoltages  in  the  electrodes  lowering 


(a) 


(b) 


(c) 


Fig.  3.  Different  flow  fields  for  fuel  cells,  (a)  Channel  structure;  (b)  interdigitated  structure;  (c)  meander  structure. 
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Fig.  5.  Boundary  cases  for  the  flow  distribution  in  a  meander  structure  combined  with  a  diffusion  layer.  Case  (A)  shows  an  ideal  flow  through  the 
meander  channel.  For  high  permeabilities  the  flow  type  (B)  can  occur  where  the  flow  passes  only  the  diffusion  layer. 


the  current  density.  In  the  following,  parameters  character¬ 
izing  the  flow  homogeneity  are  defined. 

In  a  system  consisting  of  a  meander  structure  and  a  diffu¬ 
sion  layer  there  are  two  borderline  cases.  The  first  assumes 
low  permeabilities  of  the  diffusion  layer.  Then  the  overall 
flow  is  guided  through  the  meander  only.  The  second  bor¬ 
derline  case  describes  a  high  permeability  of  the  diffusion 
layer  compared  to  the  meander.  This  results  in  the  total  flow 
going  through  the  diffusion  layer.  The  two  boundary  cases 
are  shown  in  Fig.  5. 

The  mean  specific  volume  flow  i>Spec  is  equivalent  to  the 
inlet  volume  flow  V  divided  by  the  cell  area  Acen.  This  pa¬ 
rameter  will  be  used  for  the  calculation  of  the  flow  homo¬ 
geneity  by  comparing  it  with  the  local  flow  rates. 

By  integrating  the  local  volume  flow  F1““nder  (1/h)  over 
the  entire  length  of  the  meander  channel  Lm,  the  specific 
volume  flow  in  the  meander  channel  i>spec  is  obtained  as 

7,- meander  1  1  /  .'/meander  j  t  i  . 

%ec  =7 - 7—/  ^loeal  dLM  (2.1) 

Acell  Lm  Jo 


to  the  overall  flow  as  defined  in  Eq.  (2.3)  can  be  used  in  the 
same  time  as  an  expression  for  the  flow  homogeneity  i,//m  : 

t  meander 

ifM  =  XM  =  V-^—  (2.5) 

r’spec 

If  the  flow  is  passing  only  in  the  meander,  the  local  flow 
VC-"der  is  constant  and  according  to  Eq.  (2.1)  the  parameter 
i//m  has  the  value  1.0.  If  the  flow  only  passes  in  the  diffusion 
layer  the  flow  takes  mainly  the  shortest  connection  between 
cell  inlet  and  outlet.  Then,  the  corners  are  expected  to  be  less 
supplied  with  reactants  (Fig.  5)  whereas  the  areas  between 
inlet  and  outlet  are  well  supplied.  This  causes  the  presence 
of  insufficiently  and  sufficiently  supplied  areas  on  the  MEA. 
To  take  this  into  account,  the  maximum  value  of  Xdiff  is  set 
to  0.5  for  the  borderline  case  of  the  flow  passing  only  in  the 
diffusion  layer. 

rdiff 

if  diff  =  7^  (2.6) 

r’spcc 


In  contrast  to  the  flow  in  the  meander  (Fig.  5A),  where  the 
flow  velocity  is  always  in  the  same  order,  the  flow  velocity  in 
the  diffusion  layer  strongly  varies  across  the  electrode  area 
(Fig.  5B).  In  the  region  of  the  inlet  and  the  outlet,  the  flow 
cross-section  is  substantially  smaller  than,  e.g.  in  the  middle 
of  the  cell.  The  total  specific  flow  through  the  diffusion  layer 

-rdiff  •  ,  •  ,  , 

i>spec  is  obtained  by  mass  conservation: 

rdiff  r  r meander 

^spec  ^spec  ^spec  (2-2) 


The  fractions  xm  and  Xdiff  are  reflecting  the  ratio  of  the 
average  flow  in  both  the  meander  and  diffusion  layer  to  the 
total  flow. 


XM  = 


r  meander 
V 

‘'spec 

Fspec 


rdiff 

V 

.  v spec 

and  Xdiff  =  = - 

Fspec 


(2.3) 


Due  to  mass  conservation,  the  fraction  Xdiff  of  the  flow 
passing  the  diffusion  layer  is  thus: 


Xdiff  =  1  —  XM 


(2.4) 


In  the  following  we  define  the  flow  homogeneity  for  both 
the  meander  and  diffusion  layer.  The  ratio  of  meander  flow 


By  using  xm  and  Xdiff  as  weighting  factors,  the  overall 
flow  homogeneity  i jr  results  in 

if  =  XM^M  +  Xdiff  ^diff  (2-7) 


respectively. 


(2.8) 


This  is  a  semi-empirical  expression  valid  for  a  meander 
flow  field  combined  with  a  porous  diffusion  layer.  Of  course, 
if  the  bipolar  plate  itself  is  porous  like  porous  carbon  mate¬ 
rials  with  a  meander  structure  manufactured  inside  then  the 
resulting  flow  can  be  described  equivalently. 

The  borderline  cases  are  summarized  in  Table  2  for  the 
two  borderline  cases  mentioned  above.  Due  to  insufficient 
flow  supply  to  the  corners  the  overall  homogeneity  is  is  0.5 
for  the  flow  passing  only  in  the  diffusion  layer.  In  contrast, 
flow  passing  only  in  the  meander  leads  to  a  overall  homo¬ 
geneity  t/t  of  1 . 

Flow  homogeneity  helps  to  evaluate  different  flow  fields, 
but  it  is  not  the  sole  criterion  for  the  evaluation  of  flow  fields 
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Table  2 

Parameters  of  flow  homogeneity  for  borderline  cases 
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Flow  passing  only 

in  the 

Meander 

Diffusion  layer 

Tdiff  , t 

0 

i 

^specAspec 

Pdiff 
t meander 
^spec  fv  spec 

0 

1 

0.5 

0 

if m 

1 

0 

yf  (overall  homogeneity) 

1 

0.5 

or  for  a  decision  whether  a  flow  field  is  suitable  or  not. 
Another  important  parameter,  especially  for  the  cathode  of 
fuel  cells,  is  the  pressure  drop  which  requires  compressor 
work  reducing  the  total  system  power.  In  general,  a  flow 
field  is  a  compromise  between  flow  homogeneity  and  low 
pressure  drop.  In  our  simulation  studies,  the  pressure  drop 
will  be  taken  into  account  too. 

2.2.  Reynold’s  numbers  of  the  anode  and  the  cathode 

An  important  issue  for  the  flow  distribution  in  the  mean¬ 
der/diffusion  layer  compound  is  the  ratio  of  the  flow  con¬ 
ductivity  with  respect  to  the  drag  coefficient  of  the  meander 
channel  and  the  permeability  of  the  diffusion  layer. 

Therefore,  in  a  first  step,  the  pressure  drop  characteris¬ 
tics  of  the  meander  and  the  diffusion  layers  are  determined 
separately.  To  avoid  too  high  pressure  drops,  the  meander 
structure  should  be  operated  in  the  laminar  flow  regime.  The 
Reynold’s  number  Re  is  the  decisive  flow  parameter  with: 

^  _  vdhyd  _  vdhydP  ^  ^ 

v  1 7 

with  v  being  the  flow  velocity  and  c/|lyd  being  the  hydraulic 
diameter.  The  viscosity  can  be  expressed  as  a  kinematic 
viscosity  v,  which  is  the  dynamic  viscosity  p  divided  by 
the  density  p.  Re  <  2300  indicates  the  laminar  flow.  The 
Re-number  can  be  calculated  for  the  anode  of  the  DMFC 
which  is  operated  with  a  liquid  methanol/water  mixture  and 
for  the  gaseous  cathode  of  both  DMFC  and  PEM  fuel  cells. 

The  stoichiometric  flow  rates  A.MeOH  and  Aan-  define  the 
ratio  of  the  reactants  fed  to  the  fuel  cell  to  the  amount  of 
reactants  electrochemically  converted.  To  prevent  depletion 
of  methanol  at  the  anode  or  air  at  the  cathode,  the  fuel  cell 
is  fed  with  flow  rates  X  >  1 . 

Although  the  progresses  in  DMFC  development  during 
the  last  years  lead  to  better  performance,  the  power  density 
is  significantly  lower  than  in  PEM  cells  [9].  For  a  first  es¬ 
timation  of  the  Re-number,  the  following  assumptions  are 
made: 

DMFC-anode: 

•  maximum  current  density  i  is  0.4  A/cm2; 

•  methanol  concentration  CMeOH  is  1  M; 

•  the  stoichiometric  flow  rate  of  the  methanol/water; 

•  mixture  kMeOH  is  in  the  range  of  3-8; 


•  the  cell  area  Aceii  is  100cm2. 

PEM-catliode: 

•  maximum  current  density  i  is  0.8  A/cm2; 

•  the  stoichiometric  air  flow  rate  /,lur  is  in  the  range  of 
1.5-3; 

•  the  cell  area  Aceii  is  100cm2. 


Both  flow  types  are  regarded  as  single  phase  flows.  This 
is  a  simplification  as  the  liquid  flow  of  the  DMFC  anode  is 
influenced  by  the  CCh-bubble  formation  especially  at  high 
current  densities.  At  the  cathode,  the  influence  of  the  liquid 
phase  is  relatively  small  because  the  specific  volume  of  the 
liquid  phase  is  much  smaller  than  the  gas  phase  leaving  the 
geometry  of  the  gas  channels  unchanged.  The  velocity  v  of 
the  air  in  the  meander  channel  can  be  calculated  from  the 
mass  flow  using  Faraday’s  law  and  the  geometric  parameters 
of  the  meander  channel.  When  using  N  parallel  meanders 
with  a  cross-sectional  area  Am.  the  total  volume  flow  is 
distributed  homogeneously  on  the  N  meander  channels. 


V 


AmN 


with 


M-m 


I  1 

T  _  A~P  ^ai1' 

4  F  X()2  Pair 


V  = 


1 


6  F  CMeOH 


Z-MeOH 


for  the  cathode  and 


for  the  DMFC  anode 


1000/ 


(2.10) 


In  the  equation  for  the  anodic  volume  flow  the  term  I/6F 
is  the  consumption  of  methanol  according  to  Faraday’s  law. 
The  anode  of  the  DMFC  is  fed  with  a  methanol/water  so¬ 
lution  which  has  a  concentration  CMeOH  expressed  in  mol/1. 
Typical  concentrations  are  in  the  range  from  0.2  to  1  M 
depending  on  the  current  density,  the  operating  temperature 
and  membrane  properties  as  thickness,  methanol  diffusivity 
and  other  parameters  [10,11].  To  prevent  methanol  deple¬ 
tion  in  the  channels,  the  methanol/water  mixture  is  fed  with 
a  stoichiometric  flow  rate  /,\icoii  >  1.  In  comparison  to 
the  cathode,  this  flow  rate  can  be  significantly  higher  as  the 
energy  demand  for  the  methanol/water  pump  is  generally 
lower  than  that  for  the  air  compressor.  A  typical  range  for 
^MeOH  is  4-10. 

Geometric  parameters  of  the  meander  channel  are  the 
height  Iim  and  the  width  /?m-  The  cross-sectional  area  Am  of 
each  meander  channel  is 

Am  =  hyibyi  (2.11) 


The  hydraulic  diameter  c/hyd  of  the  meander  channel  is 
calculated  by  the  cross-sectional  area  Am  and  its  circum- 
fence  Um- 


dhyd  — 


4Am 

Um 


4Z?m/?m 
2(/?m  +  Iim) 


(2.12) 


By  combining  the  Eqs.  (2.10)  and  (2.12),  Re  can  be  ex¬ 
pressed  for  both  the  anode  the  cathode  with  N  representing 
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Table  3 

/te-number  for  different  operating  conditions 


Type 

Temperature 

in  (°C) 

25 

50 

100 

Liquid  DMFC-anode 

(^-MeOH  =  3-8,  I  =  40  A) 

232-618 

534-1425 

703-1874 

Gaseous  PEM-cathode 
(Xair  =  1.5-3,  /  =  80  A) 

989-1979 

935-1871 

834-1667 

Dimensions  of  the  meander  channel:  =  1  mm,  h m  =  1  mm.  Methanol 

concentration:  1 M.  Active  cell  area:  100  cm2.  One  meander  channel 

(N=  1). 


the  number  of  meander  channels: 


.  1  1 

Re  =  2V - 


1 


N  v  Z?m  +  h  m 


(2.13) 


The  resulting  /fa-numbers  for  the  DMFC  anode  and  the 
PEM/DMFC  cathode  are  summarized  in  Table  3.  It  can  be 
seen  that  the  /?e-number  is  in  all  cases  smaller  than  2300.  At 
high  current  densities  and  higher  air  flows,  the  f?<?-numbers 
are  in  the  range  1600-2000.  An  additional  increase  in  air 
flow  could  turn  the  flow  into  turbulence.  This  is  an  indication 
that  at  high  current  densities  and  high  air  flow  rates  the  total 
flow  should  be  divided  on  at  least  two  channels  to  avoid 
turbulent  flow. 


2.2.1.  Pressure  drop  in  the  meander 

The  fte-number  calculated  above  is  an  important  parame¬ 
ter  for  calculating  the  pressure  drop  in  the  meander.  In  gen¬ 
eral  the  meander  structure  consists  of  straight  sections  which 
are  connected  in  the  corners  where  the  flow  is  redirected. 
The  pressure  drop  in  the  straight  section  of  the  meander  is 

Ap„  =  r— V  (2.14) 

^hyd  ^ 


with  r  being  the  drag  coefficient  for  laminar  flow: 


64  64/u 

r  =  cp —  =  tp - 

Re  pdhydV 


(2.15) 


The  correction  factor  tp  fits  the  drag  coefficient  r  to 
the  geometry  of  the  specific  hydraulic  diameter.  For  the 
quadratic  cross-section  of  the  channels  used  in  the  experi¬ 
ments  < p  is  0.88  [12].  In  addition  to  the  pressure  drop  in  the 
meander  channel,  the  redirection  of  the  flow  in  the  corners 
of  the  meander  causes  further  pressure  drops.  They  can  be 
calculated  by 

Apredirect  0“  (2.16) 


with  C  being  in  the  range  of  1-2  depending  on  the  geometry 
of  the  corners  (square-edged  or  rounded). 


3.  Experimental 

3.1.  Permeability  measurements  of  the  diffusion  layer 

To  determine  the  permeability  k  of  the  diffusion  layers  the 
experimental  set-up  shown  in  Fig.  6  was  used.  The  material 
sample  is  mounted  between  two  endplates  made  of  perspex. 
One  endplate  is  provided  with  concentric  inlet  and  outlet 
channels.  The  volume  flow  V  measured  by  a  mass  flow 
meter  is  pressed  through  these  concentric  channels  and  the 
pressure  drop  A p  is  detected.  A  defined  thickness  .v  of  the 
diffusion  layer  can  be  adjusted  by  distance  spacers  which 
can  be  mounted  additionally  between  the  endplates  which 
servers  to  investigate  different  compression  stages  of  the 
diffusion  layers. 

The  experimental  set-up  is  suitable  for  measurements  with 
air  and  with  water.  When  using  water,  the  air  mass  flow 


inner 

channel 


outer 
diameter  R, 


flow  direction 


out 


inner 

diameter  R;, 


Fig.  6.  Experimental  set-up  for  the  determination  of  the  pressure  drop  of  diffusion  layers. 
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controller  is  replaced  by  a  pump  and  the  volume  flow  is 
measured  by  weighing  the  outgoing  water  stream. 

The  flow  velocity  passing  the  porous  diffusion  layer  can 
be  calculated  using  Darcy’s  law: 

fdiff  =  — V/7  (3.1) 

>1 


The  measured  volume  flow  V  is 
V  = 


kdp 

-Inr - i 

i]  dr 


By  integration  of  Eq.  (3.2) 
y  rRoi 
2ns  „ 


rR out  dr  k  CPoa' 

'JRm  r  il  J  Pm  P 


(3.2) 


(3.3) 


the  permeability  coefficient  k  is  obtained  with  /%,  —  pout 
being  the  pressure  drop  A p: 


i±ln(y 

A p  2ns  V  ^  in  / 


(3.4) 


The  permeability  A:  is  a  geometric  parameter  expressed  in 
the  unit  m2. 

The  diffusion  layers  to  be  measured  were  manufactured 
in  house.  To  obtain  high  power  densities  an  optimized  mass 
transport  to  the  catalyst  region  is  necessary.  The  prepara¬ 
tion  of  the  diffusion  layers  and  the  influence  of  porosity 
on  the  power  density  is  one  of  the  most  important  topics 
and  has  been  investigated  in  detail  [13-18].  In  our  experi¬ 
ments  both  the  DMFC-anode  and  PEM-cathode  gas  diffu¬ 
sion  layers  were  prepared  by  mixing  carbon  powder  with 
finely-dispersed  PTFE.  With  regard  to  the  different  mass 
transport  properties  of  the  electrodes,  the  PTFE  content  is 
in  the  range  of  15%  in  the  anode  to  provide  a  sufficient 
transport  of  liquid  methanol  and  water  to  the  catalyst  layer 
and  40%  in  the  cathode  for  an  improved  hydrophobicity  for 


Table  4 

Influence  of  the  height  s  and  the  medium  (water/air)  on  the  permeability 
of  a  diffusion  layer 


Gap  width  s  (mm) 

Permeability  k  (m2) 

Air 

Water 

0.2 

1.05E-11 

- 

0.3 

1.41E— 11 

3.34E-12 

0.4 

1.62E-11 

6.60E-12 

0.5 

3.33E-11 

2.34E-11 

At  the  same  time  the  height  s  indicates  the  compression  state.  Structural 
data  are  the  same  as  in  Fig.  7. 


product  water  removal  [16].  The  mixture  was  then  applied 
to  carbon  cloth  with  loading  of  ca.  5  mg/cm2  with  a  subse¬ 
quent  sintering  process  at  350  °C.  An  important  difference 
between  anode  and  cathode  diffusion  layers  is  the  pretreat¬ 
ment  of  the  carbon  cloth:  in  the  case  of  anode  diffusion 
layers  the  carbon  cloth  is  untreated,  whereas  for  cathode 
diffusion  layers  the  carbon  cloth  is  hydrophobized  by  PTFE 
with  a  loading  of  ca.  1  mg/cm2.  Fig.  7  shows  the  measured 
pressure  drop  as  a  function  of  the  volume  flow  for  different 
gaps  widths  s.  The  pressure  drop  is  nearly  proportional  to 
the  flow  rate  except  for  the  case  of  the  widest  gap  (0.5  mm) 
which  is  in  the  same  order  as  the  thickness  of  the  diffusion 
layer.  This  can  be  explained  by  the  formation  of  bypass 
flows  between  endplates  and  diffusion  layer  if  the  endplates 
of  the  measuring  device  (Fig.  6)  are  not  sufficiently  pressed 
on  the  diffusion  layer.  Consequently,  at  smaller  gap  widths 
with  a  better  contact,  this  effect  cannot  be  further  observed. 

Although  the  permeability  is  a  geometric  property  of  the 
diffusion  layer,  it  is  different  for  water  and  for  air  as  shown  in 
Table  4.  Due  to  the  surface  tension  of  water,  the  pores  of  the 
diffusion  layer  become  filled,  changing  the  pore  geometry 
and  decreasing  the  permeability.  During  the  operation  of  a 
fuel  cell,  the  diffusion  layers  are  partially  wetted,  therefore 


Fig.  7.  Pressure  drop  with  variation  of  the  compression  of  the  diffusion  layer.  In  uncompressed  state  (5  =  0.5  mm)  the  pressure  drop  is  ca.  eight  times 
smaller  than  in  the  fully  compressed  state  (s  =  0.2  mm).  Composition  of  the  diffusion  layer:  30%  PTFE,  70%  XC72  on  hydrophobized  carbon  cloth. 
Measurements  made  in  air. 
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Table  5 

Influence  of  the  PTFE-content  on  the  permeability  of  the  diffusion  layers 
PTFE  content  (%)  Permeability  k  (m2)  Ratio 


Air  Water 


15  (CC  untreated) 

1.36E-11 

2.86E-12 

4.76 

20  (CC  untreated) 

1.68E— 11 

6.22E-12 

2.70 

30  (CC  teflonized) 

1 .4  IE —  1 1 

3.34E-12 

4.22 

40  (CC  teflonized) 

3.73E-11 

9.48E-12 

3.93 

The  gap  width  is  0.3  mm  in  all  cases.  CC:  carbon  cloth. 


the  real  permeability  coefficient  is  expected  to  be  somewhere 
between  the  two  borderline  cases. 

The  PTFE  content  influences  both  the  permeabilities  for 
air  and  water  as  shown  in  Table  5.  Roughly  speaking,  an 
increasing  PTFE  content  leads  to  a  higher  permeability. 
Regarding  the  water,  this  is  an  unexpected  effect  as  one 
would  assume  lower  permeabilities  due  to  the  increased 
hydrophobicity.  The  ratio  (kajr/kwater)  is  always  in  the  range 
between  2.5  and  5  and  does  not  show  any  dependency  on 
the  PTFE  content  of  the  carbon  powder/PTFE  mixture.  This 
phenomenon  may  be  explained  by  the  slightly  different 
preparation  of  the  diffusion  layers  and  has  not  further  been 
examined. 

3.2.  Flow  visualization 

For  the  flow  visualization,  we  used  the  test  equipment 
shown  in  Fig.  8  which  is  similar  to  that  described  in  [19,20]. 
It  consists  mainly  of  a  pump  and  a  flow  field  made  of  per¬ 
spex.  The  flow  field  is  designed  in  the  same  way  as  the 
plates  of  our  laboratory  fuel  cells.  The  size  of  the  flow  field 
is  100  mm  x  100  mm.  The  pressure  loss  through  the  cell  is 
monitored  by  means  of  two  simple  manometers.  To  ensure 
a  better  visualization,  diffusion  layers  with  defined  geome¬ 
tries  made  of  polyethylene  sheets  (Fig.  9)  were  used.  In  two 
polyethylene  sheets  with  a  thickness  of  80  |xm,  channels  with 
a  defined  geometry  are  cut  (Fig.  9a).  The  resulting  geometry 
of  the  two  sheets  is  a  grid  structure  as  shown  in  Fig.  9b.  For 
a  channel  width  of  0.5  mm  and  a  channel  distance  of  5  mm 
we  measured  a  permeability  k  of  5.06E—  11m2.  This  is  in  the 


Fig.  9.  Example  of  model  diffusion  layers  made  of  polyethylene  sheets 
with  defined  geometry  (thickness  80  p,m,  channel  width  0.5  mm,  channel 
distance  5  mm).  Two  sheets  (a)  positioned  above  each  other  as  shown  in 
(b)  form  a  defined  grid  structure  with  a  permeability  k  of  5.06E— 11m2. 
This  is  in  the  range  of  real  diffusion  layers  (Table  4). 

same  order  of  magnitude  as  for  the  real  diffusion  layers  (see 
Table  4).  Furthermore,  the  usage  of  model  diffusion  layers 
extends  the  permeability  range  to  be  tested.  The  pump  of  the 
test  rig  can  be  switched  between  pure  water  and  colored  wa¬ 
ter.  A  video  camera  serves  to  record  the  resulting  flow  front. 

A  comparison  of  the  flow  distribution  for  two  different 
model  diffusion  layers  is  shown  in  Fig.  10.  The  thick  lines 
show  the  flow  front — determined  each  2  s — after  analysis  of 
the  video  sequences.  The  flow  distribution  shown  in  the  left 
of  Fig.  10  is  homogeneous,  the  whole  active  area  is  covered 
by  the  flow  which  takes  approximately  16  s.  Only  in  the  left 
and  right  side  of  the  flow  field  the  flow  has  a  slight  tendency 
towards  the  cell  outlet  which  is  not  very  important  as  this 
effect  is  limited  only  to  a  small  part  of  the  overall  active  area. 
In  contrast,  the  higher  permeability  of  the  diffusion  layer 
shown  in  the  right  part  of  Fig.  10  leads  to  an  inhomogeneous 
flow  distribution.  In  comparison  to  the  left  picture,  it  takes 
the  flow  ca.  50%  longer  (24  s)  to  cover  the  whole  area.  There 
are  two  preferred  directions  of  the  flow:  the  first  one  is  from 
the  flow  inlet  following  the  left  and  the  upper  border  of 
the  active  area,  whereas  the  second  preferred  direction  is 
the  path  following  the  bottom  and  the  right  border.  As  a 
consequence,  the  central  areas  of  the  flow  field  are  depleted. 
This  is  an  important  observation  as  consequently  the  local 
stoichiometric  flow  rate  in  the  central  area  is  much  lower 


manometer 


Fig.  8.  Test  equipment  for  flow  visualization. 
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Fig.  10.  Flow  visualization  for  different  combinations  of  meander  flow  fields  and  diffusion  layers.  The  curves  in  the  flow  field  indicate  the  flow  propagation 
velocity.  Meander  channel:  2mm  x  2mm,  active  area:  100cm-.  Permeabilities  of  the  diffusion  layers:  2.75E— 11  nr  (left)  and  3.32E—  9m2  (right).  Flow 
rate  in  each  case  1.441/h  water  with  traces  of  red  ink.  Temperature:  25  °C. 


causing  lower  local  current  densities  in  real  fuel  cells.  As 
an  estimation  based  on  the  video  sequences,  the  flow  in  the 
central  area  is  only  the  half  of  the  flow  in  the  border  areas. 
For  example,  after  9  s  approximately  the  half  of  the  whole 
cell  area  is  covered  for  the  low  permeability  of  the  diffusion 
layer  (2.75E— 1 1  m2,  left  part  of  Fig.  10),  whereas  after  the 
same  time  the  high  permeability  lowers  the  coverage  to  the 
area  after  the  same  time  by  approximately  50%  to  a  amount 
of  approximately  only  a  quarter. 

In  the  following  chapter,  we  investigate  in  more  detail  the 
influence  of  design  parameters  as  flow  field  size,  meander 
geometry,  number  of  meanders  and  permeability  of  the  dif¬ 
fusion  layers  on  the  flow  homogeneity. 

3.3.  Cell  testing  of  laboratory  scale  DMFC 

In  addition  to  the  cells  described  above,  DMFC  laboratory 
cells  using  two  different  flow  fields  shown  in  Fig.  1 1  were 
tested.  The  meander  flow  field  shown  in  the  left  consists  of 


channels  and  ribs  with  a  width  of  3  mm  each  and  a  height  of 
0.8  mm.  An  additional  channel  structure  with  a  thickness  of 
0.3  mm  is  positioned  above  the  meander  flow  field  to  provide 
reactant  flow  also  towards  the  contact  area  between  mem¬ 
brane  electrode  assembly  (MEA)  and  flow  field.  For  flow 
visualization  purposes,  this  metallic  channel  structure  is  re¬ 
placed  by  polyethylene  sheets  with  a  suitable  thickness.  On 
the  right,  the  flow  distribution  in  a  grid  structure  is  visual¬ 
ized.  The  grid  structure  is  also  provided  with  an  additional 
channel  structure.  It  can  be  clearly  seen  that  the  combina¬ 
tion  of  the  meander  stmcture  with  the  additional  channel 
structure — which  is  equivalent  to  an  additional  porosity — 
causes  an  inhomogeneous  flow  especially  in  the  middle  of 
the  cell.  The  main  flow  is  concentrated  on  both  sides  of  the 
flow  field  leaving  the  upper  middle  area  of  the  cell  unsup¬ 
plied.  The  active  area  of  each  test  cell  is  25  cm2. 

The  flow  in  the  grid  structure  is  much  more  homogeneous. 
On  the  photograph  it  can  be  seen  that  the  flow  covers  the 
whole  area  of  the  cell. 


Fig.  11.  Comparison  of  DMFC  flow  fields.  (Left)  Meander  flow  field;  (right)  grid  flow  field  with  internal  manifolds.  Both  flow  structures  are  combined 
with  a  thin  (0.3  mm)  channel  structure  which  is  positioned  directly  on  the  meander,  respectively,  the  grid. 
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3.4.  Manufacturing  of  membrane  electrode  assemblies 
(MEAs) 

The  MEAs  used  for  the  experiments  were  fabricated 
in-house.  The  membrane  electrode  assembly  (MEA)  for  our 
DMFC  test  cells  consists  of  a  polymer  membrane  (Nation 
117)  with  anode  and  cathode  catalyst  layers  and  carbon  gas 
diffusion  layers.  In  the  following,  the  main  manufacturing 
steps  of  the  MEAs  are  briefly  described.  Details  are  given 
elsewhere  [21,22]. 

Both  the  anode  and  cathode  gas  diffusion  layers  were  pre¬ 
pared  by  mixing  carbon  powder  with  finely-dispersed  PTFE 
as  described  in  detail  above.  With  regard  to  the  different 
mass  transport  properties  of  the  electrodes,  the  PTFE  con¬ 
tent  is  10%  in  the  anode  to  ensure  a  sufficient  transport  of 
liquid  methanol  and  water  to  the  catalyst  layer  and  40%  in 
the  cathode  for  an  improved  hydrophobicity  for  product  wa¬ 
ter  removal. 

The  anodic  catalyst  loading  is  4  mg/cm2  carbon-supported 
Pt/Ru  with  an  atomic  ratio  of  1:1;  cathodic  loading  is 
4  mg/cm2  Pt-black.  The  cathode  catalyst  layer  was  prepared 
by  a  spraying  method,  whereas  the  anode  was  prepared  by 
a  decal  method,  which  is  based  on  three  steps: 

1.  mixing  the  catalyst  with  a  Nafion  solution  with  subse¬ 
quent  application  on  a  PTFE  foil; 

2.  drying  process  at  80  °C; 

3.  hot-pressing  of  the  catalyst  layer  onto  the  membrane  at 
130  °C,  0.5  kN/cm2,  5  min; 

4.  removal  of  the  PTFE-foil. 

In  contrast  to  the  anode  catalyst  layer  the  cathode  was 
prepared  by  the  application  of  the  catalyst  on  the  diffusion 
layer  with  subsequent  hot-pressing  onto  the  membrane. 


Fig.  12.  DMFC  single  cell  with  an  active  area  of  25crrr.  The  end  plates 
are  made  of  stainless  steel. 

3.5.  Flow  distribution  versus  power  density 

The  MEAs  were  assembled  into  a  test  cell  shown  in 
Fig.  12.  For  the  evaluation  of  the  flow  fields,  the  volt¬ 
age/current  density  characteristics  were  investigated.  Air 
was  used  as  oxidant  at  ambient  pressure  at  constant  1  Emin. 
The  methanol  concentration  was  set  to  1  M  at  a  volume  flow 
of  540ml/h.  By  means  of  a  preheater,  the  methanol/water 
mixture  (anode  inlet)  was  set  to  a  defined  temperature.  In 
addition,  the  anode  outlet  temperature  was  measured.  The 
two  temperatures  can  be  different  due  to  heat  generation 
in  the  cell  and  due  to  water  evaporation  in  the  cathode.  At 
anode  temperatures  above  60  °C,  the  anode  outlet  tempera¬ 
ture  is  significantly  lower  than  the  anode  inlet  temperature 


Fig.  13.  Comparison  of  cell  voltage/current  density  characteristics  of  laboratory  scale  DMFC  single  cells  with  different  flow  fields.  Anode  feed  temperature: 
80  °C,  volume  flow  540ml/h.  Anode  outlet  temperature:  54  °C  (meander);  64  °C  (grid  structure).  Methanol  concentration  1  M.  Air  flow  rate  1 1/min  at 
ambient  pressure. 
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Fig.  14.  Influence  of  the  temperature  on  the  power  density  for  different  flow  fields.  Anode  volume  flow:  540ml/h;  methanol  concentration:  1  M;  cathode 
volume  flow:  1 1/h  air  at  ambient  pressure. 


as  the  water  evaporation  in  the  cathode  chamber  dominates 
the  heat  balance. 

Fig.  13  shows  the  cell  voltage/current  density  characteris¬ 
tics  for  the  two  different  flow  fields  shown  above.  Compared 
to  the  grid  structure,  the  power  density  obtained  with  the 
combined  meander/channel  structure  is  in  the  order  of  only 
40%.  Due  to  the  higher  power  density,  the  temperature  in 
the  grid  structure  is  approximately  10  °C  higher  than  in  the 
meander/channel  structure  even  if  both  cells  are  provided 
with  the  same  inlet  temperature.  For  a  better  comparison 
of  the  two  cell  flow  fields,  the  power  density  is  shown  as  a 
function  of  the  cathode  outlet  temperature  in  Fig.  14.  The 
gap  in  the  two  power  densities  grows  with  rising  tempera¬ 
ture.  This  is  an  indication  for  mass  transport  problems.  An 
increase  in  the  power  density  is  accompanied  by  an  increas¬ 
ing  reactant  demand  in  the  whole  cell  area  which  can  not  be 
fulfilled  by  the  meander/channel  flow  field  causing  diffusion 
overvoltages. 


A  B  C  D  E 


4.  Simulation  study 

The  interaction  of  the  meander  and  the  diffusion  layer 
was  analyzed  by  simulation  studies.  The  program  is  written 
in  a  Fortran  95  code.  The  simulation  model  is  based  on  the 
finite  volume  element  method.  In  the  laminar  flow  regime, 
the  meander  structure  can  be  regarded  as  a  series  of  discrete 
resistances,  whereas  the  diffusion  layer  is  an  additional  re¬ 
sistance  matrix  (Fig.  15). 

The  simulation  is  based  on  the  discretization  of  the  in 
Chapter  2  defined  mass  transport  equations  similar  to  the 
method  described  in  [19].  For  each  volume  element,  the 
equations  for  the  mass  conservation  are  solved,  i.e.  the  sum 
of  the  incoming  and  outgoing  volume  flows  must  be  zero: 

V] rl  o ui  =  0  for  each  volume  element  (4.1) 

The  volume  elements  are  connected  by  the  channels  and 
the  diffusion  layers.  The  incoming  and  outgoing  volume 
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Fig.  15.  Model  of  flow  resistance  in  meander  flowfield  (left)  and  diffusion  layer  (right). 
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flows  depend  on  the  local  pressure  gradient  and  the  local 
conductivity.  For  the  laminar  flow  regime,  the  resulting  local 
pressure  gradient  is  proportional  to  the  local  volume  flow 
which  leads  to  a  system  of  linear  equations  for  the  mass 
conservation  which  can  be  expressed  as: 

Lp  =  C  (4.2) 


The  coefficient  matrix  L  contains  the  coupling  conductiv¬ 
ities  of  each  point  with  the  adjacent  volume  elements  of  the 
discretized  model  region.  The  overall  conductivity  between 
two  points  of  the  grid  is  a  superposition  of  the  conductivity 
of  the  meander  and  the  diffusion  layer  and  is  calculated 
using  Eq.  (2.14)  for  the  meander  and  Eq.  (3.1)  for  the 


\\r  =  0.48 


y/  =  0.5 


Fig.  16.  Flow  distribution  for  different  permeabilities  of  the  diffusion  layer.  The  flow  distribution  is  calculated  for  a  cell  area  of  200  mm  x  200  mm  and 
a  1mm  x  1mm  meander.  Permeabilities:  (I)  2.2E— 14m2;  (II)  4.6E— 12m2;  (III)  4.6E— 10m2;  (IV)  4.6E— 4m2.  Thickness  of  the  diffusion  layer  in  all 
cases:  0.4  mm. 
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Fig.  17.  Influence  of  the  permeability  of  the  diffusion  layer  on  the  flow  homogeneity  for  different  cell  geometries.  Results  based  on  simulation.  The 
dots  represent  the  correlation  to  the  experimentally  tested  flow  field  geometries  (see  Fig.  10).  (A)  Cell  area  200  mm  x  200  mm,  meander  cross-section 
1  mm  x  1  mm.  (B)  Cell  area  100  mm  x  100  mm,  meander  cross-section  1  mm  x  1  mm.  (C)  Cell  area  200  mm  x  200  mm,  meander  cross-section  2  mm  x  2  mm. 
(D)  Cell  area  100  mm  x  100  mm,  meander  cross-section  2  mm  x  2  mm. 


diffusion  layer.  In  a  first  step  the  vector  p  containing  the 
local  pressures  is  calculated  by  solving  (4.2).  C  comprises 
the  local  sources  and  sinks  as,  e.g.  the  channel  in-  and  outlets. 
Once  the  vector  p  is  calculated,  the  local  volume  flow  can  be 
determined  in  a  second  step  by  multiplying  the  local  pressure 
gradient  by  the  local  conductivity.  Typically,  a  system  of 
25  x  25  elements  has  been  used  for  calculation. 

4.1.  Simulation  results 

Results  of  the  simulation  studies  are  presented  in  Fig.  16. 
The  inlet  volume  flow  is  taken  as  a  reference  and  set  to  1 .0. 
Across  the  cell  area  the  relative  volume  flow  varies  accord¬ 
ing  to  the  permeability  of  the  diffusion  layer  resulting  in  a 
varying  overall  flow  homogeneity  i jr  as  defined  in  Eq.  (2.7). 
Fig.  16,  it  shows  the  orientation  of  the  meander  channel  and 
the  flow  inlet  and  outlet.  As  expected,  low  permeabilities  of 
the  diffusion  layer  lead  to  an  overall  flow  mainly  through 
flow  paths  in  the  diffusion  layer. 

Fig.  17  summarizes  the  influence  of  the  permeability  of 
the  diffusion  layer  for  four  different  cell  geometries  using 
the  simulation.  The  flow  homogeneity  i fr  is  calculated  ac¬ 
cording  to  Eq.  (2.8).  The  thickness  of  the  diffusion  layer 
was  assumed  to  be  0.4  mm  in  all  cases.  The  effects  due 
to  the  compression  of  the  diffusion  layer,  e.g.  by  cell  as¬ 
sembling  or  swelling  of  the  membrane,  are  not  taken  into 
account  for  simplification  purposes.  Consequently,  the  per¬ 
meabilities  assumed  in  the  simulation  studies  can  be  2-3 
times  higher  compared  to  real  fuel  cells  (Table  4),  i.e.  in  real 


systems  the  curves  shown  in  Fig.  17  are  slightly  shifted  to 
the  right. 

4.2.  Correlation  to  flow  experiments 

The  flow  experiments  performed  in  the  experimental 
section  can  be  correlated  to  the  simulation.  The  geometry 
of  the  test  cell  with  a  2  mm  x  2  mm  meander  described 
in  Fig.  10  complies  with  the  simulation  results  in  Fig.  17, 
curve  D,  for  both  permeabilities  of  the  diffusion  layer.  Us¬ 
ing  the  lower  permeability  of  2.75E— 11m2,  the  resulting 
flow  is  mainly  in  the  meander.  This  is  an  observation  for 
both  the  experiment  and  the  simulation.  As  stated  above, 
the  flow  has  only  in  the  left  and  right  side  of  the  flow  field 
a  tendency  directly  to  the  cell  outlet.  For  this  case  the  cal¬ 
culated  homogeneity  is  in  the  range  of  0.85.  The  higher 
permeability  of  3.32E— 9  affects  the  flow  homogeneity  sig¬ 
nificantly.  Again,  both  simulation  and  experiment  show  less 
flow  homogeneity  which  is  estimated  to  be  less  than  0.5  for 
the  experimental  flow  visualization.  The  estimation  is  based 
on  the  time  necessary  to  cover  selected  parts  of  the  flow 
field  as  well  as  on  the  estimation  of  local  flow  velocities. 


5.  Conclusion 

The  flow  distribution  in  fuel  cells  has  an  important  in¬ 
fluence  on  the  power  density.  Due  to  additional  flow  paths 
offered  by  the  diffusion  layer,  even  a  meander  structure  can 
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inhomogeneously  distribute  the  reactants  across  the  flow 
field  area.  The  flow  homogeneity  can  be  calculated  in  princi¬ 
ple  by  integrating  the  local  flow  velocities  across  the  whole 
cell  area.  Borderline  cases  are:  (a)  flow  passing  only  in  the 
meander  channel,  (b)  flow  passing  only  in  the  diffusion  layer. 
The  latter  case  is  accompanied  by  a  worse  flow  homogeneity 
approximately  a  half  of  that  of  the  former.  To  avoid  depletion 
of  reactants  in  specific  regions  of  the  cell,  the  geometry  of  the 
meander  should  be  chosen  with  regard  to  the  permeability 
of  the  diffusion  layers.  Unsuitable  combinations  of  meander 
geometry  and  diffusion  layer  properties  lead  to  even  lower 
flow  homogeneities  compared  to  the  flow  passing  only  in  the 
diffusion  layer.  High  permeabilities  require  meander  struc¬ 
tures  with  low  pressure  losses,  otherwise  the  flow  homo¬ 
geneity  decreases.  The  effects  obtained  by  simulation  studies 
have  been  proven  by  flow  experiments  with  water  and  ink. 
Experiments  using  DMFC  laboratory  cells  showed  a  signif¬ 
icant  influence  of  the  flow  distribution  on  the  power  density. 
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